Secondary protons from ultra high energy cosmic ray nuclei: is the 
Greisen-Zatsepin-Kuzmin cutoff unavoidable? 
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We discuss production of ultra high energy secondary protons by cosmic ray primary nuclei 
propagating in intergalactic space through Cosmic Microwave Background (CMB) and infrared 
(IR) radiations. Under assumption that only primary nuclei with fixed atomic mass number Ao 
are accelerated, the spectrum of secondary protons is calculated. It is found that for all Ao the 
diffuse flux of secondary protons starts to dominate over that of primary nuclei at energy E ~ 
(1 — 2) x 10 19 eV, and thus the standard Greisen-Zatsepin-Kuzmin (GZK) cutoff is produced. 

PACS numbers: 98.70 Sa, 13.85.Tp 



Ultra High Energy Cosmic Ray (UHECR) puzzle is 
one of the oldest puzzles in physics, which exists for 
more than 40 years and most probably now is close to 
being resolved. It appeared in 1966 together with pre- 
diction of the Greisen-Zatsepin-Kuzmin (GZK) cutoff 
As physical phenomenon the GZK cutoff is explained by 
production of pions by UHE extragalactic protons inter- 
acting with CMB photons. In the diffuse spectrum the 
GZK feature (steepening of the spectrum) starts at en- 
ergy (3 - 5) x 10 19 eV. The UHECR puzzle was born 
simultaneously with prediction of the GZK feature, be- 
cause already in 1966 there were known at least three 
events with energies above the GZK cutoff. With time 
the number of these events increased and in 2004 there 
were known at least 16 events with energies higher than 
1 x 10 20 eV, much above the beginning of GZK cutoff. 
Now there are strong indications that steepening of the 
spectrum compatible with GZK cutoff is observed by two 
detectors, HiRes and Pierre Auger Observatory (PAO). 
The observation of the GZK feature in differential spec- 
trum measured by HiRes @ is confirmed by the measured 
value of Ei m = 5.3 x 10 19 eV in the integral spectrum. 
E 1 / 2 is a model-independent characteristic of the GZK 
cutoff in the integral spectrum [H, whose predicted 
value precisely coincides with value measured by HiRes. 
The Auger (PAO) differential spectrum is also consistent 
with the GZK steepening 

However, below the GZK cutoff the Auger data show 
the mixed mass composition, much heavier than pure 
proton composition [5]. If primary particles accelerated 
at the sources are heavy nuclei, why GZK cutoff, which 
is a signature of protons, is observed? This problem is 
further strengthened by the Auger observation of corre- 
lation with AGN at energy E > 6 x 10 19 eV Q, which 
also requires a proton-dominated composition. 

In this paper we demonstrate that the pure nucleus 
primary component accelerated in sources results in the 
diffuse spectrum which is strongly dominated by protons 
at energy above (2 — 3) x 10 19 eV, i.e. at energy just below 
the beginning of the GZK cutoff. This proton component 
is produced by photo-disintegration of primary nuclei on 
CMB. 

We follow the paper 7] in analytic calculations of spec- 



tra. We consider the expanding universe homogeneously 
filled by the sources of accelerated primary nuclei Aq with 
maximum energy E^ x = Zq x 10 21 eV, where Aq and 
Zq are atomic mass number and charge number, respec- 
tively. The generation rate per unit comoving volume 
Qao(T,z) is given by 



Qa (T,z) 
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where z is redshift, T is Lorentz factor of accelerated nu- 
cleus, 7 g is the generation index, tun is a mass of nucleon 
and Cq is emissivity, i.e. energy generated per unit co- 
moving volume and per unit time at z = 0. In Eq. ([1]) 
we assume r m in ~ 1- 

Propagating from a source a primary nucleus Aq di- 
minishes its Lorentz factor and atomic mass number due 
to interaction with background radiations, most notably 
with CMB and infrared (IR). The Lorentz factor changes 
due to adiabatic (dT/dt = —b a d(zj) and e + e~ (dT/dt = 
— 6p air (r, z)) energy losses and atomic number due to 
photo-disintegration losses (dA/dt = — bdi S (A, T, z)). An 
important characteristic of propagation is the mean time 
of nucleus photo-disintegration ta — l/£>dis = (dA/dt) -1 . 

We study the nucleus evolution in the backward time, 
whose role is played by redshift z. We consider as initial 
state the Lorentz factor T at zq — and calculate T(z) 
for the fixed value A using the equation 



dT 

dz 



[b a d(z) 



(2) 



where dt/dz — — 1/[(1 + z)H(z)] and the Hubble param- 
eter at redshift z is H(z) — Hq^J(1 + z) 3 £l m + 57 a with 
Q m and Qa being relative cosmological matter density 
and vacuum energy density, respectively. The numerical 
solution of Eq. ||2J| gives the evolution trajectory 



T g (z) = G A (T,z ,z), 



(3) 



where the first two arguments describe the initial condi- 
tions. 

The space density of species a can be found from ki- 
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FIG. 1: Flux of primary nuclei Ao interacting only with CMB (black dotted curves labelled by "Ao") and with interaction with 
IR included (red solid curves labelled by Ao) in comparison with secondary proton flux (green dashed curves labeled by "p") for 
different values of Ao and for j g = 2.3. In case of CMB+IR the secondary protons start to dominate at energy (1 — 2) x 10 19 eV. 



netic equation 

§- t n a (T,t) - ± [b a (T,t)n a (T,t)} + ^| = Q a {T,t). 

(4) 

The solution of Eq. (01 for primary nuclei (a = Aq) ac- 
celerated at the source reads 
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where T g is given by the evolution function ([3]) , the gen- 
eration function - by Eq. JT]) and analytic expression for 
dTg/dT is given in p}- The quantity r\ takes into account 
the photo-disintegration of the propagating nucleus: 



scheme in which the secondary protons from all interme- 
diate nuclei with A < Aq are taken into account. Another 
(approximate) method developed there is based on as- 
sumption of instantaneous decay (photo-disintegration) 
of primary nuclei to Aq nucleons. Numerical computa- 
tions have shown that fluxes calculated by both methods 
coincide with very good accuracy. We use here the more 
simple method of instantaneous decay. 

Every primary nucleus emitted by a source is consid- 
ered as Aq protons and the solution of Eq. ((3]) for a = p 
with t.„ = oo reads 
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The lower limit of integration in Eq. ([5]) reflects the as- 
sumption of homogeneous distribution of the sources. 
The upper limit is imposed by factor exp(— n) accom- 
panied by condition of maximum acceleration Lorentz 
factor: Q Ao (T g ) = 0, if T g > rf n c a c x . 

Secondary protons are produced in the photo- 
disintegration chain of the primary nucleus Aq in the 
processes 7 + A — > (A — 1) + N, where 7 is a back- 
ground photon and N is a nucleon. The total flux of sec- 
ondary protons is found in Q through the computation 



where T p g (z g ) is the Lorentz factor of a proton at gen- 
eration, i.e YP{zg) = G P (T, zq, z g ), and ratio dT p /dT is 
given in 0, 0- The term e~ v in Eq. (J7J) is absent be- 
cause t p = 00. Since Qa I/-A0 (see Eq. [T]), the factor 
Aq disappears from Eq. I[7|. and the flux of secondary 
protons is universal, i.e. the same for all Aq Q. 

The calculated diffuse spectra of primary nuclei and 
secondary protons are displayed in Fig. Q] for different 
Aq, for 7 5 = 2.3 and for the cases of CMB only and 
CMB+IR with IR/optical flux taken from |[. The tran- 
sition from primary nuclei to secondary protons occurs at 
energies (1 — 2) x 10 19 eV for different Aq. The physics 
of this transition is based on universality of secondary- 
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FIG. 2: Flux of primary iron nuclei (green dashed lines labelled by "Fe") and secondary protons (blue dashed lines labelled by 
"p") for different generation indices y g . The total spectrum (sum of primary iron and secondary protons) are shown by full red 
line and compared with the Auger surface detector spectrum. The energy of intersection of iron and secondary proton spectra 
increases with y g . 



proton spectrum. In the energy range displayed in Fig.JTJ 
this spectrum is produced by photo-disintegration of 
nuclei on CMB, and the secondary-proton flux is the 
same with good accuracy in all four panels. The photo- 
disintegration of nuclei at energy E A ~ 10 19 eV is dom- 
inated by IR/optical radiation, while the produced sec- 
ondary protons have energies E p ~ Ea/A, much below, 
in case of large A, the discussed transition. Therefore, 
increasing the IR flux results in suppression of primary- 
nucleus spectrum without changing the secondary-proton 
spectrum at energies of interest. It lowers the energy of 
transition. Fig.[2]demonstrates the dependence of transi- 
tion energy on generation index 7 g and gives comparison 
of energy spectrum with Auger observations [9(. The 
calculations are performed for four values of 7 S . One 
may observe that with increasing j g from 2.1 to 2.7 the 
energy of transition increases. This effect is easy to 
understand. Let us consider the unmodified spectrum 
of primary nuclei Aq and secondary protons in the toy 
model, where Lorentz factor is conserved, and photo- 
disintegration to Aq nucleons due to interaction with 
CMB is instantaneous. The ratio of the secondary pro- 
ton spectrum and unmodified nucleus spectrum in terms 
of Lorentz factor is n p (T)/n A nm (T) — Aq. From this ra- 
tio it is easy to obtain the ratio for equal energy fluxes 
n p {E)/n A nm {E) = A\~ 1b . This ratio is suppressed by 
large 7 9 and this suppression survives after the flux of 



primaries is reduced by photo-disintegration on IR ra- 
diation. This mechanism predicts that dependence of 
transition energy on 7 g is weaker for the light primaries 
and numerical calculations confirm it. 

Fig. [2] shows that agreement with Auger spectrum is 
much better for the flat generation spectra with 7 g — 
2.1 — 2.3 than for the steep spectra. The transition energy 
for 7 9 < 2.3 is less than 2 x 10 19 eV. 

The generation spectra with j g > 2.3 are disfavored by 
the observed spectrum, though it is allowed theoretically 
in the shock acceleration [IJJ]: the distribution of the 
sources over maximum energy of acceleration E^ x or 
over luminosities makes the effective generation spectrum 
steeper at the highest energies. 

Secondary-proton dominance depends on IR flux, 
which suppresses the primary-nucleus flux. We shall 
demonstrate here that if IR flux measured at the present 
epoch z = is correct, the energy of transition remains 
the same or shifts (in case of strong evolution of the 
sources) to the lower energies. We are interested thus 
in the minimum IR flux at z > 0. We can prove the 
following general statement, valid for any diffuse back- 
ground radiation: 

In the case of the generation rate of background radi- 
ation Q(e, z) = Ke~ a (l + z) m , valid up to z max , with 
arbitrary a, m > and z max , and under the assump- 
tion that the background photons are not absorbed, the 



4 




E (eV) E (eV) 

FIG. 3: The influence of IR/optical flux at epoch z on energy of transition. The spectra of primary nuclei are displayed for 
three IR fluxes, which are the same at z = 0, but differ at z > 0. The standard one §| shown by red solid curve labelled by 
"base", the larger flux with fast evolution []| is shown by green dashed line and labelled by "fast" and the absolute minimum 
of the flux at redshift z is shown by blue dotted line (it practically coincide with the baseline curve) . The transition energy 
is given by intersection of proton flux (black dotted line) with one of the primary-nucleus curves. For minimum IR flux the 
energy of transition is the same as for the standard IR flux. 



flux of diffuse background radiation at epoch z cannot be 
lower than J z (e) — (1 + z)~ 3 / 2 Jo(e), where Jo(e) is the 
flux measured at 2 = 0. 

In Fig. [3] we illustrate the effect of variation of the space 
density of IR photons riiR,(e, z) with redshift z. Three 
models of IR fluxes are used for calculations of the flux 
of UHE primary nuclei: the baseline model of Ref. Q, 
which is considered as the standard model, the fast evo- 
lution model of Ref. @] , which gives the larger IR flux at 
higher z, and the absolute lower limit for the IR flux at 
epoch z as given above. The case of the minimum IR flux 
is very close to the standard one and does not change our 
conclusions. It is explained by the fact that the fluxes of 
UHE nuclei at E > 1 x 10 19 eV are produced at very 
small z and thus interact practically with z = flux of 
IR radiation, when the source evolution is weak. 

In conclusion, the preferential acceleration of heavy nu- 
clei in the sources is naturally accompanied by domi- 
nance of the protons in the diffuse spectrum at energy 
higher than (1 - 2) x 10 19 eV with the standard GZK cut- 



off. We studied the extreme case of nuclei acceleration 
solely, and found the dominance of secondary protons. 
In shock acceleration some fraction of accelerated pro- 
tons is unavoidably present, and this primary component 
strengthens further the proton dominance [111 ]. Another 
effect which increases the proton dominance is diffusion. 
The diffusion coefficient is inversely proportional to the 
charge number Z and thus time of nuclei propagation 
from nearby sources becomes longer and nuclei are de- 
stroyed. Variation of IR flux with z also cannot increase 
the transition energy. The only effect which increases the 
transition energy is the large effective generation index of 
UHE nuclei j g > 2.3. We conclude that the GZK cutoff 
in the observational data is compatible with heavy mass 
composition below the GZK cutoff in all most reasonable 
cases. 
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